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Abstract 

^v^ , Bifurcations of solitary waves are classified for the generalized nonlinear Schrodinger equa- 

tions with arbitrary nonlinearities and external potentials in arbitrary spatial dimensions. Ana- 
lytical conditions are derived for three major types of solitary wave bifurcations, namely saddle- 
node bifurcations, pitchfork bifurcations and transcritical bifurcations. Shapes of power dia- 
grams near these bifurcations are also obtained. It is shown that for pitchfork and transcritical 
bifurcations, their power diagrams look differently from their familiar solution-bifurcation di- 
agrams. Numerical examples for these three types of bifurcations are given as well. Of these 
numerical examples, one shows a transcritical bifurcation, which is the first report of transcriti- 
cal bifurcations in the generalized nonlinear Schrodinger equations. Another shows a power loop 
K^ phenomenon which contains several saddle-node bifurcations, and a third example shows double 

QQ , pitchfork bifurcations. These numerical examples are in good agreement with the analytical 

T^lj- ' results. 
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2 ■ 1 Introduction 



^ 



Solitary waves are spatially localized and temporally stationary (or steadily moving) solutions of 
nonlinear wave equations. Solitary waves play an important role in the understanding of nonlinear 
^ I wave dynamics and thus have been heavily studied for a wide range of nonlinear wave models 

arising in diverse physical disciplines [H |2]. When the propagation constant of solitary waves or 
physical parameters in the nonlinear wave equations changes, bifurcations of solitary waves can 
occur. Indeed, various solitary wave bifurcations in miscellaneous nonlinear wave models have 
been reported. Examples include saddle-node bifurcations (also called fold bifurcations) [21 |3l HI 
[SI m El El [U [ini [H], pitchfork bifurcations (sometimes called symmetry-breaking bifurcations) 
[IIl[l2l[l3l[Ill[l5l[l6l[ni[I3[l9], transcritical bifurcations [7J , and so on. Most of these reports on 
bifurcations are numerical. In the few analytical studies, focus was on the quantitative prediction of 
symmetry-breaking bifurcation points in the nonlinear Schrodinger (NLS) equations with symmetric 
double- well potentials UJl [TH [I6l [HI [19] and the prediction of saddle- node and pitchfork bifurcation 
points in the NLS equations with periodic potentials [11]. A general treatment of these bifurcations 
and general analytical conditions for their occurrences are still lacking at this time. 



In this paper, we systematically classify solitary wave bifurcations in the generalized NLS equa- 
tions with arbitrary nonlinearities and external potentials in arbitrary spatial dimensions. These 
generalized NLS equations include the Gross-Pitaevskii equations in Bose-Einstein condensates 
|20j and nonlinear light-transmission equations in refractive-index-modulated optical media [H [2] 
as special cases. For this large class of wave equations, we derive sufficient analytical conditions for 
three major types of solitary wave bifurcations, namely saddle-node bifurcations, pitchfork bifur- 
cations and transcritical bifurcations. In addition, shapes of power diagrams near these bifurcation 
points are also derived. We will show that the power diagram near a saddle-node bifurcation is a 
horizontally oriented parabola; the power diagram near a pitchfork bifurcation is an extra power 
curve bifurcating out from a smooth power curve on one side of the bifurcation point (like a slanted 
letter 'y'); and the power diagram near a transcritical bifurcation comprises two smooth curves 
tangentially connected at the bifurcation point. These analytical results are followed by various 
numerical examples. One example shows a transcritical bifurcation, which is the first report of tran- 
scritical bifurcations in the generalized NLS equations to the author's best knowledge. Another 
example shows double pitchfork bifurcations combined with saddle-node bifurcations, and a third 
example shows a power loop phenomenon which contains a number of saddle-node bifurdations. 
These numerical examples of bifurcations are found to be in good agreement with our analytical 
results. 



2 Preliminaries 

We consider the generalized nonlinear Schrodinger (GNLS) equations with arbitrary forms of non- 
linearity and external potentials in any spatial dimensions. These equations can be written as 

iUt + V^U + F{\U\'^,x)U = 0, (2.1) 

where V^ = d"^ /dx\ + d'^ jdx\ -|- • • • -|- 9^/9x^ is the Laplacian in the A^-dimensional space x = 
(xi, X2, • • • , xat), and F{-,-) is a general real- valued function which includes nonlinearity as well as 
external potentials. These GNLS equations are physically important since they include the Gross- 
Pitaevskii equations in Bose-Einstein condensates [20] and nonlinear light-transmission equations in 
refractive-index-modulated optical media [H [2] as special cases. Notice that these GNLS equations 
are conservative and Hamiltonian. 

For a large class of nonlinearities and potentials, these GNLS equations admit stationary solitary 
waves 

t/(x,i) = e^^*n(x), (2.2) 

where m(x) is a real localized function in the square- integrable functional space which satisfies the 
equation 

V\ - /xM -F F(u^ x)u = 0, (2.3) 

and ;U is a real-valued propagation constant. Examples of such solitary waves can be found in 
numerous books and articles (see [Il[2] for instance). In these solitary waves, /i is a free parameter, 
and u(x) depends continuously on /x. Under certain conditions, these solitary waves undergo 
bifurcations at special values of /j,. Reported examples of bifurcations in Eq. ()2.3p include saddle- 
node bifurcations [21 El El [n] and pitchfork bifurcations [IIl[l3l[l4l[l5l[l6l[I71[l8l[l9]. Transcritical 



bifurcations in this equation have not been reported yet (even though they have been found in other 
nonhnear wave models [^). 

For later analysis, we introduce the linearization operator of Eq. (j2.3p . 

Li = V2-^ + a„[F(n2,x)n], (2.4) 

which is a self-adjoint linear Schrodinger operator. We also introduce the standard inner product 
of functions, 



{f,g)= r(x)5(x)dx. (2.5) 

In addition, we define the power of a solitary wave u(x; fi) as 

/CO 
u^(x;^)(ix. (2.6) 

-oo 

This power function not only conveniently characterizes solitary-wave families, but also plays an 
important role in the stability of these waves [2]. 

Our analysis of bifurcations starts with the basic observation that, if a bifurcation occurs at 
fj, = fiQ, by denoting the corresponding solitary wave and the linearization operator as 

tto(x) = n(x; ;Uo), Lio = Li|^=^o, u=uo, (2-7) 

then the linear operator Liq should have a discrete zero eigenvalue. This is a necessary condition 
for bifurcations, hence it can be used to determine where a bifurcation might occur. This condition 
is not sufficient though. Indeed, if the function F(|f^p,x) in (12. ip does not depend explicitly on 
a certain spatial dimension Xj, i.e., the GNLS equation ()2.ip is translation- invariant along the Xj- 
dimension, then for any solitary wave u(x;/i), LiUx = 0, i.e., Li has a discrete zero eigenvalue. 
But this zero eigenvalue of Li only corresponds to a spatial translation of n(x; ^) and does not 
imply solitary wave bifurcations. More will be said on this issue in the later text (see Remark 3 in 
Sec. ED. 

In the next section, we will derive sufficient conditions for three major types of solitary wave 
bifurcations. To simplify the analysis, we will focus on the case where this zero eigenvalue of Lio 
is simple. Hence we introduce the following assumption. 

Assumption 1 Suppose at a certain propagation constant fj, = fj,Q, Liq has a zero eigenvalue. 
Then it is assumed that this zero eigenvalue of Liq is simple and discrete. 

This assumption is satisfied for almost all one-dimensional bifurcations and many higher- 
dimensional bifurcations. The case of Lio's zero eigenvalue being multi-fold (repeated) can be 
similarly treated, and that will be done elsewhere. 

Remark 1 Due to Assumption 1, the zero eigenvalue of Liq is simple and discrete, thus 
this zero eigenvalue is not embedded inside the continuous spectrum of Liq. This means that the 
solitary wave no(x) at ^ = ;Uo is not an embedded soliton [2]. This fact allows us to construct 
solitary waves in the vicinity of /i = //q by perturbation series expansions without worrying about 
continuous- wave tails beyond all orders of the perturbation expansion [2| [2T | [22 | [23] . 



Under Assumption 1, we denote the single discrete (localized) eigenfunction of Lio at the zero 
eigenvalue as tpix), i.e., 

Lw^ = 0. (2.8) 

Since Liq is a real operator, the eigenfunction ip can be normalized to be a real function. Thus -0 
will be taken as a real function in the remainder of this article. We also denote 

G{u;x) = F{u^;x)u, Gk{x) = dtG\u=uo, k = 0,1,2,3, ... . (2.9) 

These notations will be used in the next sections. 



3 The main results 

In this section, we present sufficient analytical conditions for three major types of solitary wave 
bifurcations, namely, the saddle-node bifurcation, the pitchfork bifurcation, and the transcritical 
bifurcation. In addition, power diagrams of these solitary waves near bifurcation points will also 
be described. 

First we explain what these three bifurcations are. A saddle-node bifurcation is where on one 
side of the bifurcation point /io, there are no solitary wave solutions; but on the other side of hq, 
there are two distinct solitary wave branches. These two branches merge with each other as ^ — )• fiQ. 
This bifurcation is also called a fold bifurcation in the literature (following a similar practice in 
dynamical systems [24]). Examples of this bifurcation in the GNLS equation (j2.3p can be found in 
[21 |8l m [11] . A pitchfork bifurcation is where on one side of the bifurcation point /Uq, there is a single 
solitary wave branch; but on the other side of fiQ, there are three distinct solitary wave branches. 
One of these three branches is a smooth continuation of the single solution branch from the other 
side of Ho, but the other two branches are new and they bifurcate out at the bifurcation point ^o- 
As /i — > /xo, these two new solution branches merge with the smooth branch. Examples of pitchfork 
bifurcations reported so far are all symmetry-breaking bifurcations [TTl [131 [IS 113 HSl 113 HSl IT9] . 
where a smooth branch of symmetric or antisymmetric solitary waves exists on both sides of the 
bifurcation point, but two new branches of asymmetric solutions appear on only one side of the 
bifurcation point. A transcritical bifurcation is where there are two smooth branches of solitary 
waves which exist on both sides of the bifurcation point ^iq, and these solutions on both branches 
approach each other as // — )• /Uq- So far, no examples of transcritical bifurcations of solitary waves 
have been reported in the GNLS equation ()2.3p yet (to the author's best knowledge). But these 
transcritical bifurcations do exist in Eq. (12. 3p . and one such example will be presented in Sec. [5] of 
this article. 

The main result of this article is the following theorem which gives sufficient analytical conditions 
for the above three major types of solitary wave bifurcations. 

Theorem 1 Under Assumption 1, the following three statements hold. 

1. If (uo, V") 7^ and (G2, V'^) ¥" 0) then a saddle-node bifurcation occurs at ^ = ^o- When these 
two non-zero quantities have the saine (opposite) sign, the solutions bifurcate to the right 
(left) side of /i = /iQ. 



2. If {uo,i;) = {G2,^^) = 0, (1 - GsLro'^cV'') / 0, and (Gg,^'') - ^G2^^L'^^{G2^^)) + 0, 
then a pitchfork bifurcation occurs at fi = fiQ. When these two non-zero quantities have the 
same (opposite) sign, the new solution branches bifurcate to the right (left) side of /i = /io- 

3. If (no,V) =0, (G2,V^)/0, and 

(l-G2L-ino,,V'')' > (G2,^3^(G2(L-o%)2-2L-%,V'), 
then a transcritical bifurcation occurs at // = /xq- 

It is noted that under the conditions of cases 2 and 3 in this theorem, real quantities L^q uq 
and L^Q {G2ip'^), which appear in these conditions, exist (see Lemma 1 in the next section). 

Theorem 1 shows that in the generic case of {uo,ip) ^ and (G2,V'^) 7^ 0, a saddle-node 
bifurcation occurs. Pitchfork and transcritical bifurcations would arise only in more restrictive 
situations. For instance, pitchfork bifurcations generally occur only in symmetric potentials, see 
im [T3| rHj ri5| [T61 ITTI flSl flO] and Remark 2 below. Transcritical bifurcations are more rare, which 
explains why they have not been found in Eq. (|2.3p before. The above situation closely resembles 
that in finite-dimensional dynamical systems [25j. More will be said on this in the end of Sec. O 

Remark 2 An important (dominant) class of pitchfork bifurcations is the symmetry-breaking 
bifurcation. Suppose the potential in Eq. (j2.ip is symmetric, i.e., 

F(n2;-x) = F(n2;x). (3.1) 

In addition, suppose the solitary wave tio(x) has certain symmetry (even or odd in x), and the 
eigenfunction V'(x) has the opposite symmetry of uo(x) (odd or even), i.e., 

no(-x) = ±uo(x), V(-x) = tV'(x). (3.2) 

From the notation (12. 9p . we see that 

G2 = [6ud^2Fiu^;K)+4u'dl,F{u^;^)]^^^^ , 

which has the same symmetry as uq{x.). Then obviously, 

{uo,i^) = {G2,i''')=0, 

thus the conditions of Case 2 in Theorem 1 are generically satisfied. Consequently, a pitchfork 
bifurcation occurs at /x = /xq- In this case, the two bifurcated solutions ?i^(x;/x) are simply related 
as 

n~(x;/x) = n"'"(— x; /u). (3-3) 

In addition, these bifurcated solutions break the symmetry of the original uo{x) solution and are 
asymmetric, as can be seen from their asymptotic solution formulae in Eq. (I4.46P later. This 
explains why this pitchfork bifurcation is often called symmetry-breaking bifurcation in the litera- 
ture. To the author's knowledge, all pitchfork bifurcations reported so far are symmetry-breaking 
bifurcations. 

Remark 3 Suppose Eq. (12. ip is translation invariant along a certain spatial dimension Xj, i.e., 
F(|C/p,x) in (|2.ip does not depend explicitly on Xj. If this equation admits a solitary wave no(x) 



at /i = /io, then by differentiating Eq. (j2.3p with respect to Xj, we find that LiqUq^xj = 0, thus zero 
is a discrete eigenvalue of Lio with eigenfunction ip = uq^x ■ ■ In this case, simple calculations show 
that 

(no, V) = (G2, V'') = (1 - GsLro'no, V') = 0, 

and 

Thus this case does not fall into any of the three cases in Theorem 1, hence no solitary wave 
bifurcation can be predicted. This is not surprising, since a zero eigenvalue induced by translation 
invariance does not create solitary wave bifurcations in general. 

Power diagrams are important not only for displaying solitary wave bifurcations but also for 
predicting stability properties of these solitary waves |2J. The power diagrams near these three 
types of bifurcations are given in the following theorem. 

Theorem 2 Suppose Assumption 1 holds. Denoting the power of the solitary wave at the 
bifurcation point as Pq = {uo,uo), then 

1. near the saddle- node bifurcation in Case 1 of Theorem 1, power functions of the two solution 
branches u (x; ^) are 

P±(^) = Po ± Pi • (m - Mo)'/' + 0(/U - mo), (3.4) 

where the constant Pi is given by 



2. near the pitchfork bifurcation in Case 2 of Theorem 1, the power function for the smooth 
solution branch u^{x] fi) is 

PO(^) = Po + P? ■ (m - mo) + Oiii^L - fio?}, (3.6) 

where the constant P^ is given by 

P? = 2{uo,L^o^uo); (3.7) 

power functions for the two bifurcated solution branches n (x;/i) are 

P±(m) = Po + Pi-{fi- fio) + 0{{fi - mo)'/'}, (3.8) 

where the constant Pi is given by 

Pi = 2(.o,Lro'no) + Hl-G,L^,\o,r)- 

^ '° ^ (Gg,^'^)- 3(02^^2^ Lro'(G2V2))' ^ ^ 



3. near the transcritical bifurcation in Case 3 of Theorem 1, power functions for the two solution 
branches are 

P±(/i) = Po + A • (/" - Mo) + Pt ■ (^ - f^o? + 0{(/U - Mo)'}, (3.10) 

where the constants Pi and Pg ^^^ given by 

Pi = 2{uo,L-^\o), Pt = 2{uo,uf) + {uf,uf), (3.11) 

with ti]^ specified in Eq. ()4.58p . Ug being particular solutions to Eq. (|4.59p . and 5i in 
(Ii38]) - (ig]l given in (lOT]) . 

This theorem shows that the power diagram near a saddle-node bifurcation is a horizontally 
oriented parabola. Near a pitchfork bifurcation, power curves of the three solution branches are 
all linear functions of /u. In addition, the two bifurcated solution branches u (x; /x) have the 
same power slope at the bifurcation point. In fact, in the dominant case of symmetry-breaking 
bifurcations discussed in Remark 2, power curves P (//) of the two solution branches u (x; /i) 
are identical for all ^ both near and not near the bifurcation point, i.e., P^(/i) = P~(/i), due to 
the relation ()3.3p . It is also important to notice that the smooth solution branch n''(x;/x) and 
the bifurcated solution branches u^{'x.;fi) have different power slopes at the bifurcation point, i.e., 
P^ ^ Pi, because the numerator in the second term of Eq. (13. 9p is nonzero for pitchfork bifurcations 
(see Theorem 1). Thus, the power diagram near a pitchfork bifurcation looks like a slanted letter 'y'. 
Near a transcritical bifurcation, power slopes of the two solution branches at the bifurcation point 
are the same, but their curvatures are different in the generic case. Thus the power diagram near 
a transcritical bifurcation comprises two smooth curves tangentially connected at the bifurcation 
point. These features of the power diagrams (for pitchfork and transcritical bifurcations) differ 
significantly from their familiar solution-bifurcation diagrams, and these differences are illustrated 
schematically in Fig. 1. 

The upper row of this figure plots the deviation values u{xq; fi) — uo(xo), as a function of /i, 
between solitary waves m(x; /i) away from the bifurcation point and the solitary wave uo{x) at the 
bifurcation point at a representative (fixed) xq position. These curves are drawn using the leading- 
order perturbation series solutions (14.12|) (for the saddle-node bifurcation), (|4.23p and (J4.46P (for 
the pitchfork bifurcation), and (j4.68p (for the transcritical bifurcation), which we will derive in the 
next section. Notice that these deviation diagrams closely resemble the corresponding bifurcation 
diagrams (of the same names) in finite-dimensional dynamical systems ^25j. The lower row of Fig. 
1 plots the associated power diagrams for the bifurcations in the upper row. These power curves are 
drawn using the power-function's asymptotic formulae (j3.4p . (j3.6p . (j3.8p and (j3.10p in Theorem 2. 
Notice that the power diagram of the pitchfork bifurcation has a double-branching structure rather 
than the familiar triple-branching structure, and the power diagram of the transcritical bifurcation 
has a tangential-intersection structure rather than the familiar 'x'-like crossing structure. These 
power-diagram behaviors have no counterparts in finite-dimensional dynamical systems, and they 
should be borne in mind when identifying solitary wave bifurcations in the GNLS equations (j2.ip . 



4 Proofs of the main results 

To prove the main results in Theorems 1 and 2, the following lemma will be used. 
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Figure 1: Schematic plots for solitary wave bifurcations (upper row) and the associated power 
diagrams (lower row). Column (a): saddle-node bifurcation; column (b): pitchfork bifurcation; 
column (c): transcritical bifurcation. The upper row shows the deviations u(xo; ^) — uq{xq) versus 
/i at a representative xq position. These plots are drawn using the perturbation-series solution 
(j4.12p for (a), (j4.23p and (j4.46p for (b), and (|4.68p for (c). The power diagrams in the lower row 
are drawn using the asymptotic power-function formula (13. 4p for (a), (13. 6p and (13. Sp for (b), and 
()3.10p for (c). Blue and red colors in columns (b,c) represent different solution branches. 



Lemma 1 Suppose /(x) is a localized function. Then under Assumption 1, the linear 
inhomogeneous equation 

Lio^ = f (4.1) 

admits localized solutions (p if and only if the inhomogeneous term / is orthogonal to the homoge- 
neous solution V', i-e., 

(^,/)=0. (4.2) 



This lemma is a direct consequence of the Fredholm Alternative Theorem. It can also be 
proved by expanding the localized function /(x) and the solution (j){x) into the complete set of 
eigenfunctions of the Schrodinger operator Liq and then solving for </>(x) directly. 



In the later text, the orthogonality condition ()4.2p will be called the solvability condition of the 
inhomogeneous equation (14. ip (for the existence of localized solutions). 

Proof of Theorem 1 We will use the constructive method to prove this theorem. Specifically, 
we will explicitly construct solitary wave solutions, in the form of perturbation series expansions, 
which exist near /i = /ig under the conditions of this theorem. It will be shown that these pertur- 
bation series solutions can be constructed to all orders. The existence of available solitary wave 
solutions near fi = fiQ will readily reveal the type of bifurcations at fi = fiQ. 

Case 1: saddle-node bifurcations 

Here we consider the first case of Theorem 1, and show that under its conditions {uo^ip) / 



and (G2, "0^) 7^ 0, there exist two solitary- wave branches on only one side of fi = fiQ, which merge 
with each other as /i — t- ^o- We will also show that no other solitary wave solutions can be found 
near // = //q- Hence a saddle-node bifurcation occurs here. 

The solitary waves which exist near /x = /xq in this case have the following perturbation series 
expansions 



u{^;i^) = "^{f^ - fJ'o)''^'^Uk{:s.) 



fc=0 

= uo(x) + (^- /Uo)^/^ni(x) -h (/x-/io)^i2(x) H . (4.3) 

Inserting this expansion into Eq. (j2.3p and recalling the notations ()2.9p . we get the following 
sequence of equations for Uk at order [fi — ^o)'^ > k = 0,1,2,3, ... : 

V'^uo-l^oUo + F{ul,x)uo = 0, (4.4) 

^10^1 = 0, (4.5) 

L10U2 = no - -^G2ul, (4.6) 

L10U3 = ui- G2U1U2 - ^Gsuf, (4.7) 

Lion4 = U2- —^G2 (^2 + 2nin3) - ^Gauf U2 - -r\'^^u\, (4.8) 

The equation (j4.4p for uq is satisfied automatically since uq is a solitary wave of Eq. ()2.3p at /x = /xq • 
The Ml solution to Eq. (I4.5p . under Assumption 1, is 

wi = ^'iV', (4.9) 



where h\ is a constant. The M2 function satisfies the linear inhomogeneous equation (14. 6p . Due to 
Lemma 1, Eq. (j4.6p admits a localized solution for U2 if and only if 

(V^,Mo-^G2n2) = 0. (4.10) 

Inserting the u\ solution (j4.9p into this orthogonality condition and recalling the assumptions of 
Case 1, we find that 

Thus, we get two h\ solutions itry which are opposite of each other. Inserting the corresponding u\ 
solutions (j4.9p into ()4.3p . we then get two perturbation series solutions of u{x; fi) as 

n±(x; /i) = no(x) ± ry(/i - /io)'/'V'(x) + 0(^ - /^o)- (4.12) 

If {uo,ip) and (G2,'0^) have the same sign, then r; is real. Recalling that lio(x) and ipi^) ^^^ both 
real localized functions, we see that these two perturbation series solutions (j4.12p give two real- 
valued (legitimate) solitary waves when ^u > fiQ, but not when // < //q- On the other hand, if {uQ,ip) 



and {G2,tp^) have the opposite sign, tj is purely imaginary. In this case, the perturbation series 
solutions ()4.12p give two real- valued solitary waves when ^ < ^q, but not when /i > /ig- 

Next we show that the two real localized perturbation series solutions ()4.12p . which exist on 
only one side of fx = ^q, can be constructed to all orders of {ji — /xq)^'^. Let us first consider the M2 
equation (j4.6p . When 61 is selected from Eq. (j4.1ip . the orthogonality condition (j4.10p is satisfied. 
Thus by Lemma 1, localized solutions for U2 exist. Since the inhomogeneous term and the linear 
operator Lio of (j4.6p are both real, these localized U2 solutions can also be made real. Let us denote 
one of such real localized U2 solutions as n2, i.e., 

U2 = L^Q (no - -rfG2tp^ 

then since ip \s a homogeneous localized solution of ()4.6p . the general localized solution of (j4.6p is 

U2=U2 + b2'(p, (4.13) 

where 62 is a constant to be determined. 

Now we proceed to the M3 equation ()4.7p . Inserting the ui and U2 solutions ()4.9p and ()4.13p 
into (j4.7p . we get 

L10U3 = hU- G2U2^ - jfiG^i,^ - b2G2i^A . (4.14) 

By Lemma 1, this equation admits localized U3 solutions if and only if its right hand side is 
orthogonal to the homogeneous solution tp. Imposition of this orthogonality condition yields the 62 
value as 

{l-G2U2-V^G-ii^V?,\,i>^) ,.^.. 

{G2,tV 

which is a real constant. Notice that with this 62 value, the solution n2(x) in (j4.13p is the same 
for both choices ztr/ of 61 in the ui solution ()4.9p . thus «2(x) is the same for both branches of 
the perturbation series solutions n^(x;^) in (j4.12p . With the 62 value (|4.15p . Eq. (j4.14p admits 
localized solutions 

U3 = bi {U3 + b^iP) , (4.16) 

where U3 is a real- valued localized solution of Eq. (j4.14p but without the 61 factor on its right hand 
side, and 63 is a constant. This 63 will be determined from the solvability (orthogonality) condition 
of the U4 equation (|4.8p and can be found to be real. Note that n3(x) and 63 are also the same for 
both branches of perturbation series solutions ti^(x; /i). 

Proceeding to higher orders and using the method of induction, we can readily show that all 
even terms U2n are of the form 

U2n = U2n + b2n1p, n = l,2,..., (4.17) 

and all odd terms are of the form 

U2n+1 = bl [U2n+1 + b2n+ltp] , n = 1, 2, . . . , (4.18) 
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where ti2n(x) and n2n+i(x) are certain real localized functions, and 62^, ^2n+i are certain unique real 
constants. In addition, n2n(x), 'U2n+i(x), and 62n+i are the same for both branches of perturbation 
series solutions m^(x; /i). Thus, by denoting U2n+i = U2n+i + ^2n+i'0) we have 

ufn = ■"2n(x), nf„_^i = ±r/52n+i(x). (4.19) 

Inserting these M2n ^^^ ^2n+i solutions into (14. 3p . we obtain two perturbation series solutions for 
u(x;/x), to all orders of (/x — fio)^''^, as 

oo 
U^{x;fJ,) = Uo(x) + ^(/i-/Xo)"M2n(x) 

n=l 

±r/(/. - ^o)'/' I V'(x) + f;(A^ - /io)"52n+i(x) I . (4.20) 

These two solutions exist on only one side of /i = /xq and are real and localized. The side of their 
existence depends on whether rj in (j4.1ip is real or imaginary. When // — t- //q, "U (x; /i) — )• no(x), 
thus ti^(x;/i) approach each other and merge at the bifurcation point. 

Lastly, we show that except the above two solitary wave branches which exist on only one side 
of the bifurcation point, we can not find other solitary wave solutions near this bifurcation point. 
For instance, if we look for smooth solitary- wave branches which exist on both sides of /x = fiQ, 
then their perturbation expansions should be 

m(x; fi) = uo(x) + (/U - ;Uo)ui(x) + (/x - /io)^n2(x) -\ . (4.21) 

When this expansion is substituted into p.3p . the 0(1) equation is still ()4.4p which is satisfied. At 
0(/x — ^o)) we get the equation for ui as 

LiQUi = uo. (4.22) 

Under conditions of Case 1, ('0, uq) ^ 0. Thus by Lemma 1, Eq. (I4.22P can not admit any localized 
solution for ui. This means that solitary waves with the perturbation expansion (j4.2ip can not 
exist in this case. We have also searched solitary waves near fi = fj,Q m other perturbation series 
expansions, and could not find such solutions either. Thus a saddle-node bifurcation occurs at 

fl = flQ. 

Case 2: pitchfork bifurcations 

Now we consider the second case of Theorem 1. We will show that under conditions of this 
case, a smooth branch of solitary waves exists on both sides of ;U = /xq- In addition, two other 
solitary- wave branches exist on only one side of /x = /xq- As /x — )• /xq, all these solution branches 
approach the same solitary wave uo{x). Thus a pitchfork bifurcation occurs at ^ = ^q. 

(i) We first construct the smooth branch of solitary waves which exists on both sides of /x = /xq. 
These solitary waves have the following perturbation series expansion 

oo 

u°(x; /x) = J^(/x - /Xo)'=Ufc(x). (4.23) 

fc=0 
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Inserting this expansion into Eq. (j2.3p , we get the following sequence of equations for ut at orders 
(m-^o)', A: = 0,l,2,3,...: 

V'^UQ-^xoUQ + F{ul,^)uQ = Q, (4.24) 

Lio^i = no, (4.25) 

^^10^2 = ui- -^G2u\, (4.26) 

1 „ 

Lions = U2- G2U1U2 - TTyGsni, (4.27) 

1 11 

Lion4 = ns - —G2 (n| + 2nin3) - —G2,u\u2 - -^^G^ul, (4.28) 

The equation (j4.24p for no is satisfied automatically. Under conditions of Case 2, (f/', no) = 0. Thus 
by Lemma 1, the solvability condition for the ni equation (|4.25p is satisfied, hence this equation 
admits localized solutions 

m = 2i + biijj, (4.29) 

where 

ni = L^o^no (4.30) 

is a real localized particular solution to Eq. (I4.25|) . and 61 is a constant to be determined. Inserting 
this ni solution into the U2 equation ()4.26p . we get 

Lions =ui- hj2ul + 6iV'(l - G2U1) - ^6?G2^^ (4.31) 

By Lemma 1, the solvability condition of this U2 equation is that its right hand side be orthogonal 

to the homogeneous solution ■0- Under conditions of Case 2, (G2,V''^) = 0. Thus this solvability 

condition gives 

,1 

^2' 

Since the inner product on the left side of this equation is nonzero under conditions of Case 2, this 

equation yields a unique 61 value as 



61(1 - GsL-ino, V') = {-G2ul - uu^l,). (4.32) 



.1 = ^^-"^/^-/^'^^ (4.33) 

(l-G2Lro^no,V2)' ^ ' 

which is a real constant. Hence a real localized solution for ni has been obtained. 

With the above 61 value, the solvability condition of the n2 equation (I4.3ip is satisfied. Thus 
this equation admits a real localized particular solution n2, and its general solution is 

n2 = n2 + &2V', (4.34) 

where 62 is another constant to be determined. 

Inserting this n2 solution into the ns equation ()4.27p . this equation becomes 

Lions = 62V'(1 - G2ni) + ^2(1 - Gsm) - ^Gsn?. (4.35) 
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By Lemma 1, the solvability condition of this u^ equation is that its right hand side be orthogonal 
to ip. Utilizing the ui solution (I4.29P and the conditions of Case 2, we see that 

(^'(l - Gsni), V) = (1 - G2L^^uo,^P^) / 0. (4.36) 

Thus the solvability condition of Eq. (14.350 yields a unique real 62 value, 

{U2i^-G2ui)-G3uy3\,ip) 

' (V^(l-G2^l),^) 

hence a real localized U2 solution ()4.34p is then obtained. At this 62 value, Eq. ()4.35p admits a real 
localized particular solution u^, and its general solution is 

U3 = U3 + bsij, (4.37) 

where 63 is another constant to be determined. 

Pursuing this calculation to higher orders, it is easy to see that for any n > 2, the Un solution 
is of the form 

Un = Un + bnip, (4.38) 

where Un is a real localized particular solution of the n„ equation, and 6„ is a constant to be 
determined from the solvability condition of the Un+i equation. The Un+i equation is always of the 
form 

LioUn+l = (1 -G2Ui)Un + Tn+l{uo,Ui, . . . ,'U„_i;x), (4.39) 

where Tn+i is some real function which depends on the already-obtained real localized solutions 
uo,ui, . . . , Un-i as well as x. Inserting the Un solution (14.380 into (14.390 and utilizing Eq. (|4.36p . 
the solvability condition of ()4.39p is met at a unique real 6„ value, hence a real localized u„ solution 
(j4.38p is obtained. Meanwhile, since the solvability condition of (j4.39p is met, a real localized 
particular solution Un+i exists, and the general Un+i solution is of the form (I4.38P with the index n 
replaced by n + 1. This process then repeats itself. Hence a real- valued and localized perturbation 
series solution (I4.23P is constructed to all orders, and it gives a branch of real-valued solitary waves 
'u*^(x; /i) which exists on both sides of /i = /io and depends smoothly on //. 

(ii) Next we construct two additional solitary wave branches which exist on only one side of 
/i = /io and merge with the above smooth solution branch at n = hq. These additional solitary 
wave branches have the following perturbation series expansion 

oo 

u(x; /i) = ^(/x - ;Uo)'=/2^fc(x). (4.40) 

fc=0 

This perturbation series is of the same form as ()4.3p in Case 1. Thus when this perturbation series 
is substituted into Eq. (j2.3p . the resulting equations for u^ are the same as (j4.4p - ()4.8p before. But 
since the conditions for Case 2 are different from those for Case 1, solutions Uk for the perturbation 
series here will differ from those in (14.31). as we will show below. 



First, the equation (j4.4p for uq is satisfied automatically since uq is a solitary wave of (|2.3p at 
/i = /io- The solution ui to Eq. (j4.5p . under Assumption 1, is 

Ml = blip, (4.41) 
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where 61 is a constant to be determined. Inserting this ui solution into the U2 equation (j4.6p . we 
get 

L10U2 = uo- \hlG2^^. (4.42) 

Due to conditions of Case 2, both uq and G2ip'^ are orthogonal to ip. Thus by Lemma 1, both 
L^Q uq and Lj~q {G2ip'^) exist and are certain real localized functions. Hence the solution U2 to Eq. 
(I02]l is 

U2 = L^o^no - ^blL^^\G2iJ^) + bsV', (4.43) 

where 62 is another constant to be determined. Inserting these ui and U2 solutions into (j4.7p . the 
1(3 equation is 

L10U3 = 61 |(1 - G2L-^'uo)iJ - ^bj [Gsi;^ - 3G2V'Lro'(G'2V'')] - &2G'2V''| • (4.44) 

In view of the conditions of Case 2, the solvability condition of this u^, equation yields the 61 value 
as 

Two 61 values ±1^ are obtained which are opposite of each other. Inserting the corresponding ui 
solutions (j4.4ip into (|4.40p . we then get two solutions u^{'x.;fi) as perturbation series 

u±(x; /u) = no(x) ± zy(/i - /Uo)^/V(x) + 0(/^i - ^0), (4.46) 

where ly is given in (j4.45p . If the numerator and denominator under the square root of ()4.45p have 
the same sign, then v is real. In this case, two real localized perturbation series solutions (j4.46p 
are obtained when fj, > fj,Q. If the numerator and denominator in (j4.45p have the opposite sign, 
then u is purely imaginary. In this case, two real localized perturbation series solutions (|4.46p are 
obtained when /j, < hq. 

Next we show that the two real localized perturbation series solutions ()4.46p . which exist on 
only one side of // = /xq, can be constructed to all orders of (/i — /io)^'^. With the choice of 61 values 
in ()4.45p . the solvability condition of the u^ equation ()4.44|) is met, thus the U3 solution is 

u-i = 61 [u3 - h2L^^{G2ilJ^) + hili\ , (4.47) 

where u^ is a real localized function which satisfies the equation 

L10S3 = (1 - G2Lro'no)^ - ^1/2 [GsV-^^ - 3G2V'iro'(G2^')] , (4.48) 



and 63 is a constant to be determined. Inserting these ui,U2 and U3 solutions into (|4.8p . the U4 
equation becomes 

LiQUA = 62 I (1 - G2L-^^UQ)i: - h^l [03^-^ - 3G2V'Lro' (G2^')] | 

-\hlG2iJ^+Fi{uo,^,hl^)-bihG2iJ^, (4.49) 
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where J-4 is a real localized function which depends on uq, ip, bf and other already-obtained real 
functions (such as u^). Utilizing the 61 formula (I4.45P as well as conditions of Case 2, the solvability 
condition of this U4 equation is met at the unique real 62 value, 

^-^^'^^ (4.50) 



2(l-G2Lro%,V2)" 

When this 62 value is inserted into (|4.43p . a real localized U2 solution is then obtained. Notice that 
this 62 is the same for both choices ±1/ of 61 in the ui solution (I4.4ip . thus n2(x) is the same for 
both branches of the perturbation series solutions n^(x;/i) in (|4.46p . 

For the 62 value given in ()4.50p . the solvability condition of the ii4 equation (j4.49p is satisfied, 
thus this equation admits the following localized solution 

U4 = U4-u%L{^{G2i^^) + b4i^, (4.51) 

where u^ is a real localized function which satisfies the u^ equation (j4.49p but without the last (63) 
term, and 64 is a new constant to be determined. 

Starting from n > 5, the Un equation can be derived from (|2.3p and the expansion (j4.4Up . and 
is all of the form 

LioU2n+l = U2n-1 - G2{uiU2n + U2U2n-l) - -Gzulu2n-1 + bi'H2n+l{uo,1p,bl,X), Tl > 2, (4.52) 
Ll0U2n+2 = U2n - G2{uiU2n+l + U2U2n) " -^G3ulu2n + ^2n+2(^i0, V', ^1, x), n > 2, (4.53) 

where Ti2n+i and T-l2n+2 are real localized functions which depend on uq, ip, 6f and other already 
fully determined real quantities (such as us, u^, 62, etc). Using the method of induction as well as 
conditions of Case 2, we can show that all Un solutions are of the form 

U2n+l=bl[u2n+l-b2nL^oiG2lp'^) + b2n+l1p] , n > I, (4.54) 

U2n+2 = U2n+2 " l^'^b2n+lL^o {G2'lp'^) + &2n+2V', n > I, (4.55) 

where U2n+i and 22n+2 are certain real localized functions, and b2n+i, b2n+2 are real constants 
which are determined uniquely from the solvability conditions of the ^271+3 and U2n+4 equations. 
We can also show that U2n+i, U2n+2, ^2n+i and 62,1+2 depend on 6f as a whole and are thus the 
same for both solution branches u (x;//). Inserting these solutions into the perturbation series 
(j4.40p . we obtain two branches of solitary waves ti^(x; /i) whose perturbation series expansions are 

00 
ii=^(x; /i) = tio(x) + ^(^ - /Uo)"'n2„(x) 

n=l 

±i.(^ - ^0)'/' I ^(x) + ^(/i - ^o)"S2n+i(x) i , (4.56) 

where real localized functions U2n are given by ()4.43p and (j4.55p . real localized functions U2n+i are 
as U2n+i in ()4.54p but without the 61 factor, and z^ is given in ()4.45p . These two real solitary waves 
exist on the side of ^ > /io {fi < no) when u is real (purely imaginary). When fj. ^ fiQ, they both 
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approach ito(x), thus these u^(x;/i) solution branches merge with the smooth n'^(x;/i) solution 
branch in (14.230 at ^ = fiQ. 

The existence of the smooth solution branch u^{x; fi) in ()4.23p on both sides of ^ = /io as well 
as two additional solution branches u (x;/i) in (|4.56p on only one side of /x = /io indicates that a 
pitchfork bifurcation occurs at /j, = hq. 

Case 3: transcritical bifurcations 

Now we consider the third case of Theorem 1. We will show that under conditions of this case, 
two smooth branches of solitary waves exist on both sides oi n = hq, and these branches intersect at 
/i = /io where solitary waves on the two branches become identical. Thus a transcritical bifurcation 
occurs at /i = /io- 

In this third case, we seek solitary wave solutions which exist on both sides of /i = /xq and 
depend smoothly on /i near /i = /iQ- The perturbation series expansion of such solutions is 

oo 

u(x; /i) = ^(/x - /io)''ufc(x). (4.57) 

fc=0 

The form of this expansion is the same as (j4.23|) in Case 2, thus the equations for u/^ are also the 
same as (I4.24p - (l4.28p before. However, the solutions to these equations will differ from the previous 
ones in Case 2 due to different conditions of the present case. 

The uq equation (j4.24p is satisfied automatically since uq is a solitary wave of Eq. ()2.3p at 
/i = /io- Under conditions of Case 3, the solvability condition of the ui equation (I4.25p . {uo,ip) = 0, 
is met. Thus by Lemma 1, localized ui solutions of the form 

ui = Lj"(/no + hi) (4.58) 

are admitted. Here L^qUq is a real and localized particular solution to Eq. ()4.25p . and 6i is a 
constant to be determined. Inserting this ui solution into the U2 equation ()4.26p . this equation 
becomes 

LioU2 = L^o^uo - ^G2{Lio^uof + 6i(l - GaLro^uo)^ - ^^fCsV''- (4.59) 

The solvability condition of this equation gives the following quadratic equation for bi : 

{G2,^P'')bl - 2(1 - G2L^^uo,iJ^)bi + {G2{L^^uof - 2L^^\o,4^) = 0. (4.60) 

Under conditions of Case 3, the coefficient of the bf term in this quadratic equation is nonzero, and 
A = (1 - GsLio^no, , V')' - {G2,i)^){G2{L^o'^of - 2L^^uo,i^) > 0. 

Thus this quadratic equation admits the following two real roots 

{G2,i)-^) 

For each of these two 6i values, a real localized ui solution (j4.58p is obtained. In addition, a real 
and localized particular solution U2 to the U2 equation (j4.59p exists, hence the U2 solution is 

U2 = U2 + b2ilJ, (4.62) 
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where 62 is a new constant to be determined. 

Inserting the above U2 solution into the ^3 equation (I4.27p . we get 

^10^3 = 62(1 - G2Ui)i^ + (1 - G2Ui)u2 - ^Gsn?. (4.63) 

UtiUzing the ui solution ()4.58p and the 61 formula (j4.6ip . we find that 

((1 - G2Ui)^lj,^lj) = (1 - GaL^o'^o.V^') - 61(^2, V'') = ^^ + 0. (4.64) 

Thus the solvability condition of Eq. ()4.63p yields a real constant 62 as 

, ((1 - G2Ux)u2 - G3^f/3!, f) 

' ((1 - G2ni)V^, V') 

For this 62 value, the solvability condition of the U3 equation (j4.63p is satisfied, thus this equaiton 
admits a real localized particular solution S3, and the general uj, solution is 

^^3 = ^3 + 63^^, (4.65) 

where 63 is another constant to be determined. 

Pursuing this calculation to higher orders, it is easy to see that for any n > 2, the ti„ solution 
is of the form 

Un = Un^ bnip, (4.66) 

where n^ is a real localized particular solution of the Un equation, and 6„ is a real constant to be 
determined from the solvability condition of the u„+i equation. The Un+i equation is always of the 
form 

LioUn+l = (1 - G2Ui)Un + J^n+l{uQ,Ui, . . . ,'U„_i;x), (4.67) 

where J-'n+i is some real function which depends on the already-obtained real localized solutions 
uq,ui, . . . , Un-i and x. Inserting the u„ solution (j4.66p into (j4.67p and in view of Eq. ()4.64p . the 
solvability condition of (I4.67P then yields a unique real value for the constant 6„. 

In the above solution process, since 61 can take either one of the two real roots b^ in (j4.6ip . 
ui in ()4.58p then has two corresponding solutions u^ . These two ui solutions cascade up to higher 
orders, and thus all Un functions have two solutions u^. Consequently, two real- valued and localized 
perturbation series solutions 



n± 



(x;/i)=uo(x) + X;(^-Mo)'n±(x) (4.68) 



k=l 



are obtained to all orders, and they provide two branches of real- valued solitary waves u (x; //) 
which exist on both sides oi fi = fiQ and depend smoothly on //. When fi — )• ;Uo, both u^(x;//) 
approach uo(x), thus these two solution branches intersect at /i = ^o- As a result, a transcritical 
bifurcation occurs at ^ = ^q. This completes the proof of Theorem 1. D 

Next, we prove Theorem 2 on power diagrams near bifurcation points. 
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Proof of Theorem 2 The power formula ()3.4p of saddle-node bifurcations can be derived 
easily from the perturbation series solutions ()4.3p and the ui solution (I4.9p with 5i given by Eq. 
()4.1ip . The power formula (j3.6p for the smooth solution branch u^(x;ij,) in a pitchfork bifurcation 
can be derived readily from the perturbation series solutions (I4.23P and the ui solution (I4.29p . To 
derive the power formula (|3.8p for the two bifurcated solution branches in a pitchfork bifurcation, 
we substitute the ui , U2 solutions in (I4.4ip and (I4.43P into the expansion (14.56P , and find that the 
power function is given by ()3.8p . where 

Pi = 2{uo,L^o^uo) + bl [(V', V') - {uo,L^^{G2i^^))] , (4.69) 

whose value is the same for both bifurcated branches. Since Lj~q is self-adjoint and Lj~q uq exists 
(by Lemma 1), this Pi coefficient can then be rewritten as (13. 9p (here the 6i formula (I4.45P is 
also used). The power formula (j3.10p for transcritical bifurcations can be derived easily from the 
perturbation series solutions (I4.68P and the ui, U2 solutions (I4.58p . (I4.62p . D 



5 Numerical examples of solitary wave bifurcations 

In this section, we present numerical examples for these three types of solitary wave bifurcations, 
and compare them with the analytical results presented in Theorems 1 and 2. So far, examples of 
saddle-node and pitchfork bifurcations of solitary waves have been reported in the GNLS equations 
(|2T]1 with various potentials and nonlinearities [HIHllHKIIKIlOliailllllZlIISKIS]. Here we will 
present some new examples of saddle-node and pitchfork bifurcations in the GNLS equations which 
exhibit interesting novel features. In addition, we will present the first example of transcritical 
bifurcation in these GNLS equations. 

Example 1: combined saddle- node and double-pitchfork bifurcations 

The first example we choose is the one-dimensional GNLS equation (j2.ip with a symmetric 
double-well potential and cubic-quintic nonlinearity: 

iUt + U,, - V{x)U + |C/|2[/ - 7|[/|4[/ = 0, (5.1) 

where the symmetric double-well potential V{x) is taken of the form 

V{x) = -Vo [sech^(x + xq) + sech^(x - xq)] , (5.2) 

Vo > is the potential depth, 2x0 is the separation between the two wells, and 7 > is the 
coefficient of the quintic nonlinearity. Notice that the cubic and quintic nonlinear terms in (|5.ip 
have the opposite sign, and the quintic term induces a self-defocusing effect which counters the 
self-focusing effect of the cubic term. One may also view this opposing cubic-quintic nonlinearity 
as a Taylor-series approximation to the saturable nonlinearity in photorefractive crystals J26]. The 
parameter values in the above GNLS model are chosen as 

Vo = 2.8, xo = 1.5, 7 = 0.25. (5.3) 

Solitary waves in Eq. ()5.ip are sought of the form ()2.2p . where u{x) is a real localized function 
satisfying the equation 

Uxx - ^J'U- V{x)u + u^ - ju^ = 0. (5.4) 

18 



When u{x) is infinitesimal, the hnear Schrodinger operator of Eq. (j5.4p admits a positive symmetric 
discrete eigenfunction at eigenvalue fi ^ 1.7896. This eigenmode is the ground state of the underly- 
ing double- well potential. From this linear (infinitesimal) eigenmode, a family of positive symmetric 
solitary waves bifurcates out. The power curve of this symmetric-soliton family is shown in Fig. 2 
(blue curve in the upper left panel). We have computed the spectra of the linearization operator 
Li for these solitary waves, and found that their spectra contain a simple zero eigenvalue at three 
locations marked by letters 'A,B,C' on the power curve. This is evidenced in the upper right panel 
of Fig. 2, where the Li-spectra of solitary waves at these three locations are displayed. Notice that 
at locations 'A,B', the second largest eigenvalue of the spectrum is zero, while at location 'C, the 
largest eigenvalue is zero. At these three locations, solitary waves uo{x) and eigenfunctions ■0(x) of 
the zero eigenvalue in Li's spectra are plotted in the lower row of Fig. 2 (as solid blue and dashed 
red curves respectively). Notice that eigenfunctions 'il){x) at points 'A,B' are anti-symmetric, while 
the eigenfunction at point 'C is symmetric. 
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Figure 2: Bifurcations of solitary waves in Example 1. Upper left: the power diagram; the blue 
curve is for the family of symmetric solitary waves, and the red curve is for the family of asymmetric 
solitary waves which bifurcate out from points 'A, B' through double pitchfork bifurcations. Upper 
right: Li's spectra for solitary waves at bifurcation points 'A,B,C' of the power diagram. Lower 
row: solitary waves uo{x) and eigenfunctions Tp{x) of Li's zero eigenvalue at bifurcation points 
'A,B,C' (the eigenfunctions are normalized to have unit amplitude). 



At these 'A,B,C' points, zero is a simple discrete eigenvalue of Li, thus Assumption 1 is met 
and Theorem 1 applies. In addition, for the present example, 



G2 = 6uo 



20jul 



G, = 6 



60 jul. 
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Now we use our analytical criterion in Theorem 1 to determine if and what bifurcations occur at 
these points. 

At points 'A,B', it is easy to see from symmetry that 

In addition, when the eigenfunction ip is normalized to have unit amplitude (see Fig. 2, lower row), 
we find numerically that at point 'A', 

(1 - G2L^o^uo,iP^) = -4.9313, (Gg,^^) - 3(G2^^L^o'(G'2V'')) = -58.4035; 

and at point 'B', 

(1 - G2L^oUo,ij'^) = 23.9913, (Ga,^'^) - 3(G2V^i^^o^(G2V^)) = -110.9244. 

Then according to Theorem 1, pitchfork bifurcations occur at both 'A' and 'B' points. In addition, 
the new (asymmetric) solitary waves bifurcate out on the right side of point 'A' and on the left side 
of point 'B'. 

At point 'C, we find that 

{uo,ip) = 6.4879, (G2,V'^) = -21.0632, 

thus according to Theorem 1, a saddle- node bifurcation occurs at this point. In addition, the 
bifurcated solutions appear on the left side of point 'C. 

These analytical predictions of bifurcations prove to be completely correct. Specifically, at points 
'A,B', symmetry-breaking pitchfork bifurcations occur. The two bifurcated asymmetric solitary 
waves n^(x; //) are related to each other by a mirror reflection in x, i.e., u'^{—x; fi) = u~{x; jx), and 
their power curves (which are identical) are displayed as the red line in Fig. 2 (upper left panel). 
Notice that these bifurcated solutions appear on the right side of point 'A' and on the left side of 
point 'B', as predicted by the above analysis. To illustrate solution profiles before and after these 
bifurcations, we focus on point 'A'. The power diagram near this bifurcation point is amplified from 
that in Fig. 2 and shown in Fig. 3 (first panel from the left). Notice that the power curves near 
this bifurcation point are linear functions of /i, in agreement with Theorem 2 and Fig. 1(b). We 
have also compared the slopes of the power curves at point 'A' in Fig. 3 with the analytical power 
slopes in Eqs. (|3.6p - (l3.9p . and found excellent agreement. At three locations 'a,b,c' on the two 
sides of the bifurcation point 'A' in the power diagram, profiles of the solitary waves are displayed 
in Fig. 3(a-c) respectively. Solutions in Fig. 3(a,b) are symmetric and lie on the symmetric branch 
of the power diagram (blue line), while the two solutions in Fig. 3(c) are asymmetric and lie on 
the asymmetric (bifurcated) branch of the power diagram (red line). Notice that on the left side of 
the bifurcation point, there is a single (symmetric) solitary wave (see Fig. 3(a)); but on the right 
side of the bifurcation point, there are three solitary waves, one symmetric (see Fig. 3(b)) and the 
other two asymmetric (see Fig. 3(c)). These behaviors of the pitchfork bifurcation agree fully with 
our analytical results as well as the schematic plots in Fig. 1(b). 

It is interesting to observe from Fig. 2 (upper left panel) that the asymmetric soliton branch 
starts out from point 'A' and terminates at point 'B', thus it appears and then disappears through 
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Figure 3: Power diagram and profiles of solitary waves near the pitchfork bifurcation point 'A' in 
Fig. 2. First panel: power diagram; (a,b,c) solitary waves at locations marked by the same letters 
on the power diagram. 

double pitchfork bifurcations. In between, its power curve exhibits a 'S' shape, indicating that 
double saddle-node bifurcations also occur on this asymmetric branch. These features of bifurca- 
tions are quite novel for the GNLS equations (j2.ip . In a different nonlinear wave system, namely 
nonlinear saturable couplers, a similar double pitchfork bifurcation also exists [12j. 

At point 'C of Fig. 2, we have found that a saddle- node bifurcation occurs as predicted. This 
is already obvious from the power diagram in Fig. 2, which shows that the power curve turns 
around at this point. The power diagram near this saddle-node bifurcation point 'C is amplified 
and shown again in Fig. 4. This numerical power curve is compared with the analytical saddle-node 
power formula (|3.4p and complete agreement is obtained. At two locations 'a,b' of the power curve 
below and above the bifurcation, profiles of the solitary waves are displayed in Fig. 4(a,b). These 
solutions are all symmetric, and their amplitudes vary while going through the bifurcation. 




Figure 4: Power diagram and profiles of solitary waves near the saddle-node bifurcation point 'C 
in Fig. 2. Left panel: power diagram; (a,b) solitary waves at locations marked by the same letters 
on the power diagram. 

Example 2: power loop phenomena 

Our second example is still the GNLS equation with opposing cubic and quintic nonlinearities, 

iUt + U^x - V{x)U + \U\^U - 0.15|C/|^C/ = 0, (5.5) 

but the double-well potential V{x) is now asymmetric instead: 

V(x) = -3.5 sech2(x + L5) - 3 sech2(x - 1.5). (5.6) 

This potential is displayed in Fig. [5]^a). As usual, solitary waves in this equation are sought of the 
form (j2.2p . where u{x) is a real localized function. We find that in this system, there exist a family 
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of positive solitary waves whose power curve forms a closed loop. This power loop is displayed in 
Fig. [5]^b). This power loop shows that this family of solitary waves has a non-zero minimal power 
and a finite maximal power, and it exists over a finite propagation-constant interval. In addition, 
four saddle-node bifurcations are clearly visible on this loop. We have checked that at these saddle- 
node bifurcation points, the bifurcation conditions in Theorem 1 (Case 1) are all satisfied. At four 
locations of the power loop, three of them ('c,e,f') being saddle-node bifurcation points and the 
remaining one ('d') slightly below a saddle-node bifurcation point, profiles of the solitary waves 
are displayed in Fig. [5l[|c-f ) . It is seen that the energy of these solitary waves resides primarily in 
the shallower (right) well of the potential. Thus this family of solitary waves is different from the 
family of ground-state solitary waves in this system, whose energy resides primarily in the deeper 
(left) well of the potential. One may notice that this power loop in Fig. [5]^b) self-crosses itself in 
the middle (above point 'd'). This power-curve crossing does not signal a transcritical bifurcation 
however, because as /i approaches this crossing point along the two intersecting curves, the solitary 
waves do not approach each other. This power loop phenomenon of solitary waves has not been 
reported before in the GNLS equations (12. ip (to the author's knowledge), but a similar momentum 
loop phenomenon for solitons sitting on constant backgrounds has been discovered in the externally 
driven NLS equation [27j. 




Figure 5: Power-loop phenomenon in Example 2 (i.e., Eq. (15.51) ). (a) The asymmetric double-well 
potential V{x) in Eq. (j5.6p : (b) the power loop; (c-f) profiles of solitary waves at locations marked 
by the same letters on the power loop of (b). 

Example 3: transcritical bifurcation 

Our last example is the GNLS equation with competing cubic, quintic and seventh-power non- 
linearities, 

iUt + U^x - V{x)U + \U\'^U - 0.2\U\^U + K\UfU = 0, (5.7) 

where V{x) is the same asymmetric double-well potential (j5.6p as in Example 2, and k is a real 
constant. In this example, a transcritical bifurcation of solitary waves is found at 



K = Kr^ 0.01247946. 



(5.^ 
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The power diagram of this bifurcation is shown in Fig. ^h). We see that two smooth solution 
branches, namely the upper C1-C2 branch and the lower di-d2 branch, tangentially connect at the 
bifurcation point {ho,Pq) « (3.28,14.35). Profiles of solitary waves at the marked ci,C2,di,d2 
locations on this power diagram are displayed in Fig. [Hc-d). Notice that these solutions are 
close to each other since the ci,C2,di,d2 locations are near the bifurcation point (/xo,-Po)- As fi 
approaches this bifurcation point, we find that these solitary waves along both the lower and upper 
power branches approach each other, confirming that this is a transcritical bifurcation. Notice that 
the power diagram in Fig. [6|^b) agrees with the analytical power formula (]3.10p of transcritical 
bifurcations (see also the schematic power diagram in Fig. 1(c)). In addition, we have checked the 
conditions of transcritical bifurcations in Theorem 1 (Case 3), and found them satisfied here as 
well. 
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Figure 6: Transcritical bifurcation in Example 3 (see Eq. (|5.7p ). (a) The asymmetric double-well 
potential V{x) in this example; (b) the power diagram; (c) profiles of solitary waves at locations ci 
(solid blue) and C2 (dashed red) of the upper power curve in (b); (d) profiles of solitary waves at 
locations di (solid blue) and d2 (dashed red) of the lower power curve in (b). 



What would happen to the bifurcations in the above three examples if the potential or the non- 
linearity in those underlying GNLS equations is slightly perturbed? We have numerically studied 
this question and found that in Example 1, when the nonlinearity or the potential is slightly and 
arbitrarily perturbed (including perturbations to make the double- well potential ()5.2p asymmetric) , 
the saddle-node bifurcations (at point 'C of the symmetric-soliton branch and two others on the 
asymmetric-soliton branch in Fig. [2]) always persist. For the two pitchfork bifurcations in this 
example (at points 'A,B' of Fig. 2), if the perturbed potential is still symmetric, then these pitch- 
fork bifurcations would survive; but if the perturbed potential becomes asymmetric, then these 
pitchfork bifurcations are destroyed. In Example 2, the four saddle-node bifurcations on the power 
loop of Fig. [5] always persist under weak perturbations in the nonlinearity or the potential. In 
Example 3, the transcritical bifurcation is extremely sensitive and is destroyed under generic small 
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perturbations to the system (such as when k ^ Kc). From these numerical results, we conclude that 
saddle-node bifurcations are generic and robust under weak perturbations to the system; pitchfork 
bifurcations are generally reliant on a symmetric potential; and transcritical bifurcations are very 
fragile and generally disappear under perturbations. These behaviors are consistent with similar 
statements on these bifurcations (based on the bifurcation conditions) below Theorem 1. 



6 Summary and discussion 

In this paper, we classified solitary wave bifurcations in the generalized NLS equations ()2.ip with 
arbitrary nonlinearities and external potentials in arbitrary spatial dimensions. Sufficient analyti- 
cal conditions were derived for three major types of solitary wave bifurcations, namely saddle-node 
bifurcations, pitchfork bifurcations and transcritical bifurcations. These conditions show that the 
generic solitary wave bifurcation is the saddle-node bifurcation; the pitchfork bifurcation generally 
requires certain symmetry conditions; and the transcritical bifurcation is rare. For these bifurca- 
tions, shapes of power diagrams near the bifurcation points were also obtained. It was shown that 
the power diagram for a pitchfork bifurcation exhibits double branching rather than the familiar 
triple branching, and the power diagram for a transcritical bifurcation features two curves tangen- 
tially touching each other rather than the familiar 'x'-crossing. Numerical examples for these three 
types of bifurcations were presented as well. These examples show novel features such as power 
loops and double pitchfork bifurcations. The example of transcritical bifurcation seems to be the 
first report of such bifurcation in the generalized NLS equations (j2.ip . 



The results in this paper are important not only for a general classification and understanding of 
solitary wave bifurcations in the generalized NLS equations ()2.ip . More importantly, the bifurcation 
conditions in Theorem 1 will be the basis for a general treatment of linear stability of solitary waves 
near bifurcation points. This stability analysis lies outside the scope of the present article and will 
be reported elsewhere. We do want to say here that the stability properties of solitary waves 
near bifurcation points in Eq. ()2.ip show some important qualitative differences from those in 
finite-dimensional dynamical systems [25j. Details will be presented in a forthcoming article. 
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